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The emission of correlated electrons from an itinerant ferromagnet following the 
impact of a polarized electron beam is analyzed in terms of irreducible tensorial 
parameters that can be measured. Under favorable conditions, specified in this work, 
these parameters are related to the spin polarization in the ferromagnet. The formal 
results are illustrated by numerical studies of the polarized electron pair emission from 
a Fe(llO) surface and a novel technique for the investigation of magnetic properties 
of ferromagnets is suggested. 
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The electronic and magnetic properties of low dimensional systems with long-range mag- 
netic order, such as ultrathin ferromagnetic films and magnetic surfaces are currently under 
intensive investigations [|1]]. This is due to the fundamental and technological importance of 
such materials. Magnetic systems with reduced symmetry can be explored by a variety of 
techniques |TJ]. Low-energy spin-polarized electron spectroscopy is particularly suitable as 
the penetration depth is on the order of few atomic layers fl2|. In this method one resolves the 
quantum states of the incoming and outgoing electrons to extract the accessible information 
on the sample under investigation. On the other hand, a promising technique emerged in 
recent years where an electron pair, resolved in energy and momentum, is detected following 
the impact of an unpolarized electron beam upon a non- magnetic sample 0-§J. As demon- 
strated successfully for a variety of materials 0-|5[], the electron pair carries, under favorable 
conditions, direct signature of the Bloch spectral function which is a central quantity as far 
as the electronic structure is concerned. However, these studies 0,f| have been performed at 
higher energies (~ 20 keV) and the role of the spin polarization has not been yet addressed. 
Very recently, however, it has been demonstrated by a pioneering experiment |7| that the 
electron-pair emission depends strongly on the spin polarization of the electron beam and 
the magnetization of sample. Thus, it seems timely to inspect theoretically the low-energy 
polarized electrons emission from ferromagnets. We conclude: a) the electrons' spectra are 
quantified fully by a set of irreducible tensorial components; b) under certain circumstances 
specified below, the electron-pair spectrum is directly related to the spin-resolved spectral 
function of the surface. 

For a theoretical formulation we consider a reaction in which two electrons are simul- 
taneously emitted from a ferromagnet with a defined magnetization direction M after the 
impact of a mono-energetic spin polarized electron beam. The spins of the electrons in the 
incoming beam and in the sample are assumed to be good quantum numbers. A corre- 
sponding experiment resolves the asymptotic wave vectors of the impinging and the two 
emitted (vacuum) electrons which we label and k' x , k' 2 , respectively. No spin analysis of 
the outgoing electrons is performed. The target surface is described by the state vectors 
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|0£,a,s 2 ,m S2 ), where e is the ground state energy, m S2 is the projection of the spin s 2 of the 
ground state along a quantization axis and a denotes all other quantum numbers. The spin 
polarization of the incoming beam with projection m Sl of the electron's spin si is charac- 
terized by the density matrix Pm sl m Sl whereas the population of the magnetic sublevels of 
the state \(p t ,a,s 2 ,m a2 ) is given by the density matrix Pm 32 m 32 - The scattering probability is 
related to (atomic units, a.u., are used throughout) 

^(k^k^ki) = C 4^(K>K> m s 1 > m & a > m >a>ki> m n) 
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M*(k[,k' 2 ,m s > i ,m s > 2 ;a,m S2 ,k 1 ,m Sl )5(E f - E,) (1) 

where Ef is the final-state total energy. The initial-state total energy E^ is E^ = E^ —e whilst 
E^ is the energy of the projectile beam and C = (2ir)' l /ki. All energies are measured with 
respect to the vacuum level. The transition amplitude A4(k[, k' 2 , m s > , m s ^; a, m S2 , ki, m Sl ) 
is given by M = (Vv^k^my ,m s , |^|0e,a,s 2 ,m S2 Vkj.sim^ ) where <Pki,*im 41 is a spinor vacuum 
state describing the incoming beam. The emitted electrons with spin projections m s i^m s ' 
are represented by the state vector iVv^k'^m , ,m ,) whereas T is the transition operator. In 
Eq.(|II) the density matrices are diagonal. This is not a restriction as they can always be 
diagonalized by a rotation in the appropriate spin space. Furthermore, we adopt M as a 
joint quantization axis for si and S2. In case si and s 2 do not have a common quantization 
axis, we apply an appropriate spin rotation of the density matrix of the incoming beam. For 
convenience we express the electrons' final state in the total spin space as 
J2sm s (SMs\ s^m'^, s^m'g^ ^ , ki / ,k 2 ';5A/ s ) where S is the total spin and Ms is its projection. 

To disentangle geometrical from dynamical properties we expand the density matrices 
in state multipoles (statistical tensors) p pq ||, 

P msi m sl = E {-) pi ~' 1 ~ m ' 1 (si-m ai ;s 1 m ai \ Pl q 1 = 0)p Pl9l=0 . (2) 

Pi=0 
2s 2 

Pm S2 m S2 ( e , "hE {-) P2 - S2 - ms * (s 2 - m S2 ; s 2 m S2 \p 2 q 2 = 0)p P2q2=0 {e, a). (3) 

P2=0 
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Substituting Eqs.(0,||) into the general expression ([I]) yields 
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The decisive point is that the sum over m Sl {m S2 ) in Eq.([5D defines the component (along 
M) of a spherical tensor of rank p\ (p 2 ) This mathematical observation yields important 
information as to the transformation behaviour of Aqo P2 : Ao 1 ~°' P2 (Aqo P2_0 ) is a scalar with 
respect to spin rotations generated by Si (s 2 ), i.e. it represents spin averaged quantities in 
the s x (s 2 ) spin space, whereas the components A p ) 1 ^~ odd,P2 (k p ) 1 Q 2 ~ odd ) can be regarded as a 
spin orientation in the s x (s 2 ) spin space (for pi — 1 it is a polar vector) and hence changes 
sign upon spin reflection, i.e. A p ) ^ odd,P2 (— m Sl ) = — AQ^ odd,P2 (m Sl ) [Aq^ 2_ocM (— m S2 ) = 
— K p ) 1 Q 2 ~ odd {m S2 )\. The tensorial components with even pi values are alignment parameters, 
i.e. they describe the deviations in the spectra from the unpolarized case. The above 
formalism is easily generalized |§ to the case of strong spin-orbit coupling and/or multi- 
electron emission. For two electrons Eq.(|j) reduces to 
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The first term of the sum in Eq. (|6]) is the pair emission rate averaged over the spin orientation 
of the incoming electron beam and the spin polarization of the sample. The second term 
describes the spin asymmetry due to the inversion of the magnetization while the incoming 
electron beam is being unpolarized . The third term is the spin asymmetry in the electron- 
pair emission from unpolarized targets when inverting the spin polarization of the electron 
beam. In absence of explicit spin interactions in the transition operator T, e.g. spin-orbit 



coupling, the parameters A q and Aq'o vanish. Finally the last term of Eq.(||), the prime 
focus of the following calculations, is related to the electron-pair emission from spin-polarized 
samples by spin polarized electrons. It is a polar vector both in the Si and S2 spin spaces, 
i.e. Aq q(— m si , m S2 ) = — Ag o(m sl , m S2 ) = Aoo(m Sl , — m S2 ). The explicit forms of Aqq and 



Aq|q are derived from Eq.@ to be 



A o!o = % E E { \MK K SM S - k 1; a, I, ^ | 2 - \M(k[, k' 2 , SM S ; k 1; a, | 
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+ |-M(k' 1; k' 2 , SM S - k l5 a, T, 1r) I 2 + |-M(k' 1; k' 2 , SM,; k 1; a, |, ^) | 2 } . (8) 

The projections of the spins of the sample state and the electron beam parallel (anitparallel) 
to the quantization axis are labeled, respectively by the arrows ft (JJ-) and f (j). In the total 
spin space Eqs.(|7|,|8|) are expressed in terms of the singlet and the triplet partial cross sections, 
X( 5=0 ) and X^ s=1 \ respectively, i.e. 
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x^i^X-M-^) -x^t&X'M-,")] (9) 

SX {s=1) (k[, k 2 ; ki; a) + X {s=0 \k[, k 2 ; k i; a)] =: 2X*°*. (10) 



X( 5 =°) and A^ 5=1 ) are determined by the matrix elements, T^(k^, k' 2 ; k 1; a), of the singlet 
(S* = 0) and triplet (S = 1) transition operators T s = (1 + (— 1) 5 P 12 )T. Here P 12 is a 
permutation operator that interchanges the two emitted electrons. Thus, one obtains the 
symmetry property T^^k^, k' 2 ; kj; a) = (— l) 5 T^(k 2 , k^; k 1; a), i.e. in situations where an 
interchange of the electrons does not modify the ionization dynamics the triplet scattering 
amplitude and hence X^ 5=1 ^ = C\T^ S=1 ^\ 2 vanishes. An example will be shown below. 
Till this point the electronic and structural properties of the sample have not been yet 
specified. For perfect clean surfaces the integral over a in Eq.@ implies summation over 
the surface Bloch wave vector k 2 y and over the surface layers. The Bloch theorem imposes a 
conservation law for the surface components of the total wave vector of the emitted electrons 



K|j~ = +k^y E3' i- e - ^he change of Kjj~ from its initial value kiy + k 2 y (before the collision) 
is restricted to a multiple of the surface reciprocal lattice vector gy. This fact can be used 
to perform the integrals over k 2 |i in Eq.(^J). Therefore, in absence of spin interactions in the 
Operator T Eq.(|6|) reduces to a summation over the surface layers, indexed by I, and over 
g h i.e. 

W^2{2X tot (k[,k / 2l k 1 ,g ll ,l) [ Poo poo(e,A {l J) + p 10 p 10 (e,A l{ J)A%\<' 1 X-M,S\\,l)]S(E f - Ei)} (11) 

where Ay = K|" — gy — kiy. The "exchange scattering asymmetry" has been defined as 

A s X( s =V(k[, k 2 ; k 1; g„, /) - X(*=°)(k;, k' 2 ; kx, gy, Q 

• 3X(^)(k' 1 ,k' 2 ;k 1 ,g||,/)+X^=o)(k' 1 ,k' 2 ;k 1 ,g||,/)- 1 ) 

To calculate the terms in Eq.([ll]) the state multipoles pio and pio are needed. These 
can be obtained by inverting the relations (|2]|J). In the standard representation, the density 
operators of the beam and the surface are linearly expanded in terms of the Pauli matrices 
a as p Sl = 1 + Pi • a and p S2 = u>o(k 2 ||, I, e)(l + P 2 • cr) where iu (k 2 ||, I, e) is the spin- 
averaged Bloch spectral function of the layer I and Pi and P 2 are the polarization vectors. 
The sample polarization is given by P 2 = [w(k 2 [|, I, e, fr) — u>(k 2 |i, I, e, |L)]/[wo(k 2 ii, I, e)]. Here 
io(k 2 ||, I, e, m S2 ) stands for the spin and layer resolved Bloch spectral function. Thus we 
obtain p oPoo = [^o(k 2 ||, I, e)]/2 and pioPio = [^o(k 2 ||, I, e)] PiP 2 /2 and Eq.(|l|) reduces to 

w oc MH l > e ) xtot l l + AKEf - Ei). (13) 
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The asymmetry function A has been introduced as 

. D £, [w(A h l, e, t) - w(A h l, e, ^)] £ g|| X tot A s 5(E f - E { ) _ W(U) ~ W(it) 
A-n — — -t-t — - — — vin , xfE , w\ TTTTT^ , ttTTTXT " ( 14 J 



Ei< w (A h l>; , e) Eg;, X«*5(E f - E % ) W(U) + W(0[) ' 

Thus, for the calculation of the tensorial parameters two major ingredients are needed: 1) 
The spin and layer-resolved spectral function of the sample which can be obtained from 
the trace of the imaginary part of the corresponding Green function and 2) the matrix 
element of the singlet and triplet transition operators. Now we calculate the terms in 
Eq.(|l3"l) for a Fe(llO) surface. The Bloch spectral functions used here are provided by two 



independent calculations: 1) The scalar relativistic full-potential linearized augmented plane- 
wave method (FPLAPW) |TT| , |T2|| and 2) the full relativistic layer Korringa-Kohn-Rostoker 
method (LKKR) F° r the calculations of the transition matrix element we approximate 

the T operator by T = U sur f+U ee (l+G~ e U sur f) where U ee is the electron-electron interaction, 
G~ e is the Green function of the electron pair and U sur f is the surface scattering potential. 
For U sur f we employ, for a given layer, a linear combination of non-overlapping muffin-tin 
potentials |T0| . 

As stated above, for certain geometries, the triplet scattering amplitude vanishes due to 
symmetry and hence A s tends to —1 (cf. Eq . (|12"D ) . Thus, if a monolayer or a bulk system is 
considered the magnetic asymmetry P 2 can be scanned by determining VFQlT) and ^(jft)- 
This yields a direct (relative) estimate of the population of the spin states in the sample. For 
multilayered systems, we have to consider the weighting factor X tot in Eq. (p~4|) . An example 
is shown in Fig.l for a Fe(llO) sample. The two electrons are detected with fixed equal 
energies in the x — z plane and at symmetric positions with respect to the z direction while 
the incident beam direction is varied in the z — y plane. The experiment, in the geometry of 
Fig.l, is invariant under a 180° rotation with respect to the z direction. This rotation can be 
regarded as an interchange of k' x by k 2 and since 7 1 ^ =1 )(k / 1 , k' 2 ; ki, a) = — T( s=1 )(k 2 , k.[; ki, a) 
the triplet scattering (X^ 5=1 ^ = C\T^ s=lS) \ 2 ) vanishes. 

The energies e in Eq. fll4|) is determined by e = E[ + E' 2 — E^ where E[ and E' 2 are the 
energies of the vacuum electrons. Thus we tune E[, E' 2 and E^ such that e coincides with 
the Fermi energy Ep. Now by varying (3 = cos" 1 z • kx we scan P 2 along the T-N direction 
in the Brillouin zone, as shown in Fig.l. Alternatively one may fix the direction Aii = kni 
and image ^(e) by varying, e.g., the incident energy. For a polarized homogeneous electron 
gas one scans (as function of energy) the relative difference between the occupied density of 
states of the majority and minority bands. 

Away from the points of high symmetry (cf. Fig.l) the scattering dynamics, as described 
by X( 5=0 ) and X^ s=1 ^ become dominant. An example is shown in Fig.2 for (3 = 0. Again at 
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the T point (k^ = — k^y) the asymmetry A, and in particular its sign, is determined solely 
by P 2 . The For highly asymmetric energy sharing the scattering exchange asymmetry A s is 
small which leads to a reduced asymmetry A, as seen in Fig. 2. 
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Fig. 1: The asymmetry, as given by Eq. (|T4"D , for the emission of two equal-energy electrons 
from a magnetized Fe(110) surface following the impact of a polarized electron beam with 
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an energy of 35 eV. The total energy of the pair is fixed to E[ + E' 2 = 25.15 eV. The two 
electrons are detected in the y-z plane at symmetric position cos -1 z • k' x = 40° = cos -1 z • 
left and right to the z axis (cf. inset) and M || x. The angle of incidence (3 = cos -1 z ■ ki 
is varied in the x — z plane, as shown by the inset. In this geometry, the triplet scattering 
vanishes and A can be related to P2- The predominant contributions to A originate from 
the first and second surface layers. 

Fig. 2: The spin asymmetry A as function of the energy sharing (E[ — E 2 )/(E' 1 + E' 2 ) for 
a fixed total energy (E[ + E' 2 ) = 21 eV. The incident electron has an energy 26 eV and a 
polarization degree of ~ 65%. The same target as in Fig.l, however, we choose (3 = and 
the two electron detectors to lay in the x — z plane. As in Fig.l, the detectors are positioned 
at cos 1 z • kj = 40° — cos ^z • kj. The theoretical results are averaged over the angular 
resolution of the detectors. The experimental data are courtesy of Ref. M. The spectral 



functions are calculated within the scalar relativistic FLAPW method 11 
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